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A simple model is exhibited in which the remnant density of charged vortons is used to provide 
candidates for explaining the observed ultra high energy cosmic rays (above 10^° eV). These vortons 
would be accelerated in active galaxies and propagated through intergalactic medium with negligible 
losses of energy. The expected number density of observable events is shown to be consistent with 
extrapolation of the observations. The spectrum is predicted to be spatially isotropic while its shape 
is that of an atomic excitation-ionisation, i.e. with a few peaks followed by a continuum; there is 
also an energy threshold below which no vorton is visible. 
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INTRODUCTION 

The problem of very high energy cosmic rays [|p is still 
an open one : most models based on reasonable astro- 
physical assumptions seem to indicate a likely maximum 
value for the energy of any kind of emitted particle at 
most of ^ 10^® eV Q and the existence of the microwave 
background makes it impossible for a proton say to prop¬ 
agate with such energy on scales much larger than a few 
tens of Mpc ; this is the celebrated Greisen-Zatsepin- 
Kuz’min (GZK) cutoff |3|. An observation dated Octo¬ 
ber 1991 by the Fly’s Eye detector showed evidence for 
a cosmic ray at an energy of (3.2 ± 0.6)10^° eV Q and a 
few others reported showers above 10^° eV |Q (and well 
above the GZK cutoff). There does not seem to date 
to be any completely satisfactory suggestion for explain¬ 
ing these observations (see however Ref. |^), and it is the 
purpose of this paper to propose such a possible explana¬ 
tion that has been overlooked, namely that cosmic rays 
may consist of (as opposed to originate in) topological 
defects (as was first proposed in Ref. [||). 

A fairly common (neither confirmed nor excluded) pre¬ 
diction of particle physics models at energies beyond that 
of the electroweak scale is the existence of topological de¬ 
fects [0,^ . Among these, only cosmic strings seem to be 
consistent with cosmological constraints (with the pos¬ 
sible exception of textures which, as non-localized ob¬ 
jects, are not considered here). In particular, strings 
have been proposed as seeds for galaxy formation as well 
as sources for the microwave background fluctuations ob¬ 
served by GOBE. This requires that the strings appear 
at the Grand Unified (GUT) scale, a scenario that may 
breakdown in the case suggested by Witten in which 
they would be endowed with superconducting properties. 
In the latter case, string loops can form equilibrium con¬ 
figurations called vortons [p^,pj , whose stability is cur¬ 
rently under investigation , that would easily over¬ 

fill the Universe Qjl^, thus provoking an unobserved 
cosmological catastrophe. It has been suggested, as a 
possible way out of this problem, that cosmic strings 


might have appeared at a much lower symmetry breaking 
scale, estimated in the range ^ 10^ — 10® GeV (or 

at least that the formation of superconducting currents 
may have been postponed until this stage). 

A second reason to consider such comparatively low 
energy cosmic strings instead of the more popular GUT 
strings is that the latter have been shown |15| to yield 
a negligible flux of cosmic ray types events as compared 
to what is currently observed, whatever the mechanism 
involved. This is mostly based on the idea that, a string 
being a topologically stable object, emission of very high 
energy particles must proceed via some removal of the 
topological stability, which for a GUT string implies ei¬ 
ther extremely highly energetic phenomena (and we’re 
back to the previous problem of reaching the required 
high energies) or relatively rare situations (cusp evapo¬ 
ration, intercommutation and loop collapse). Adding up 
both cases, the expected flux turns out to be very tiny 
(roughly 10“^® times that observed). Our aim here is to 
show that it is possible in principle to consider high en¬ 
ergy rays as cosmic string loops themselves if the scale is 
not that of GUT but rather much below that scale. This, 
as we shall see, gives a flux that is comparable with ob¬ 
servations. 

The paper is organized as follows: in Section |, we 
recall the basic facts about vortons and how one may 
expect them to interact with ordinary particles. Then, 
Section is devoted to the interaction of a vorton with 
a proton at rest in order to estimate the order of mag¬ 
nitude of the various phenomena involved in the detec¬ 
tion of a 10®® eV event, including the probability that 
a vorton interacts with the atmosphere. It is predicted 
that the typical expected spectrum form should resem¬ 
ble that of an excitation-ionisation spectrum, namely it 
should consist of lines followed by a continuum. Besides, 
since the interaction probability is found to be quite low 
(the cross-section is typical of neutrino-hadron interac¬ 
tions at the same energies [IT^), w e predict important 
horizontal showers. In Section pTl| , we present a plausi¬ 
ble acceleration mechanism for vortons, which turns out 
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to be much simpler than most of the standard acceler¬ 
ation mechanisms for protons (Fermi-like mechanisms), 
without sharing their drawbacks. We then discuss prop¬ 
agation in Section 0 where it is shown that no GZK 
cutoff is to be expected for vortons. It is then argued 
that this model implies that the most energetic events 
should be distributed isotropically. Finally we summarize 
our findings in the conclusion where we also compare the 
expected properties of these rays with the anticipated de¬ 
tection capabilities of the future Auger observatory 0 - 

I. THE PHYSICS OF VORTONS 

The objects we shall now consider are vor¬ 
tons 1^,1^,^, namely loops of cosmic strings endowed 
with superconducting properties and stabilized by a cur¬ 
rent. They arise in symmetry breaking theories above 
the electroweak scale for which the vacuum manifold is 
not simply connected. In other words, if the high en¬ 
ergy vacuum is invariant under the transformation of a 
symmetry group G and the low energy vacuum under 
that of the group H, the necessary and sufficient condi¬ 
tion for the existence of cosmic strings is that the first 
homotopy group tti of the quotient G/H be nontriv¬ 
ial, Tri{G/H) 7 ^ {0}. A typical example is the scheme 
17(1) ^ {0}, which is also the scheme at work in the 
Landau-Ginzburg model for superconductivity [|^, and 
the strings are the corresponding vortices, much stud¬ 
ied experimentally [p^ . The vortices appear as infinite 
or in the form of closed loops which decay via emission 
of gravitational radiation. The symmetry breaking is 
achieved thanks to a Higgs field d) which in turn can be 
coupled to other fields, particularly charged ones, which 
we shall generically write as E. Here we consider only the 
case where this charge carrier is a scalar particle, noting 
anyway that because of the formal equivalence between 
bosons and fermions in the two-dimensional worldsheet 
generated by the string, one expects similar behaviors for 
any kind of particle. Besides, as was shown to be suffi¬ 
cient in previous work |2^, we shall only consider the 
coupling of E with the electromagnetic field as a pertur¬ 
bative effect not modifying much the overall dynamics of 
the vorton. 

Let us here summarize the basic properties of vor¬ 
tons 0 as we expect to observe them. They can be 
characterized by essentially two integer numbers, namely 
a topological one, N say, which specifies the winding of 
the current carrier phase around the loop, and a dynami¬ 
cally conserved number, Z, related in the charge-coupled 
case (to which the present analysis is restricted) with the 
total amount of electric charge Q that it holds through 
Q = Ze, with e the electron electric charge. Note that 
both N and Z can be positive or negative, although for 
N it is just a matter of convention, while in the case 
of Z, the sign is important when there is an external 
interaction (as in the case of acceleration for instance). 


Both these numbers are conserved at a classical level and 
one expects them to be of the same order of magnitude 
Z ~ ^ 100 [Q. From these, on constructs the loop 

angular momentum as J = ZN ^ 10"'’: although a bit 
hypothetical and thus often seen as “exotic”, vortons can 
actually be considered as classical objects. 

Since it can be shown that electromagnetic self¬ 
coupling is mostly negligible , one can evaluate the to¬ 
tal mass-energy of the vorton simply in terms of its char¬ 
acteristic radius r^^, its energy per unit length U ^ rri^ 
and tension T ~ as 

=27rr^(C/ + T)~m%, (1) 

with m the energy scale at which the strings are formed 
(essentially the string forming Higgs mass). Moreover, 
we shall assume in what follows that the current carrier 
mass scale TOo- is also of order m so that we can keep 
only one energy scale. Note however that the following 
analysis can be easily generalized for two different mass 
scales since in practice it is the current scale that is rel¬ 
evant in most calculations. Then, knowing the angular 
momentum to be given by |2|] = UTr y/2Tr permits 

to calculate the characteristic vorton circumference as 

27rr^ = (27r)i/2|^^|i/2('jjj.)-i/4 _ ^ 2 ) 

with a corresponding mass 

My ~ Zm. (3) 

Moreover, the mass scale is constrained: depending on 

various assumptions about the string network evolution 
and the rate of loop formation as well as the probability 
that arbitrarily shaped loops end up in vorton states, it 
can be shown that, if the vortons are stable then 

in order to avoid a cosmological mass excess (H^ < 1), 
the condition 

TO < 10® GeV, (4) 

must be satisfied. 

Having discussed the basic properties of vortons, let 
us now turn to a rough evaluation of the typical mass 
scale expected for to if those vortons were to be seen as 
cosmic rays, developing air showers. As we shall see in 
the following section, interaction between a vorton and 
whatever other particle occurs mainly via inelastic scat¬ 
tering resulting in the extraction of a trapped E particle. 
In other words, the current flowing along the string loop 
can be seen as a bunch of bound states which can be ex¬ 
cited provided the interaction energy is large enough. We 
therefore conclude at the existence of an energy threshold 
above which the typical expected spectrum should change 
qualitatively. Besides, and we now come to the second 
firm prediction of the model, since those particles form 
bound states, we expect them to show up in the form of a 
line spectrum, bound states being always quantized; this 
will be shown on a specific vorton model in section ||. 
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For the time being, what really matters is the existence 
of bound state energies the order of magnitude of which 
we shall now attempt to evaluate. 

Let AE be the variation of energy between two energy 
levels in the vorton, calculated in its rest frame, and 7 its 
Lorentz factor in the rest frame of the particle it interacts 
with (recall we are at the end interested in air showers 
occurring in the atmosphere where the particles interact¬ 
ing with the vorton, namely essentially quarks composing 
protons and neutrons, are supposed to be at rest). The 
energy e at which the interaction is then seen is obtained 
by transforming back to the particle’s rest frame, 

e ^ 'yAE. (5) 

Denoting by to the particle’s mass and requiring the in¬ 
teraction to actually take place gives 

7 TO ^ AE, ( 6 ) 

so that altogether, Eqs. (I) and d) imply 


from which the characteristic cosmic string scale can be 
deduced in the following way. From the fact that the 
angular momentum of the vorton is J ~ it is seen 
that the density of states scales like Z~‘^. Therefore, 
AE ^ to/Z^, so that using ®) and ( 0 ) yields 

TO ^ Z^V efh. ( 8 ) 

We shall now use these relations together with the ob¬ 
servations that have been realized on cosmic rays to nor¬ 
malize the energy levels.. 

Let us apply this evaluation to ultra high energy cos¬ 
mic rays so that the characteristic observed energy 
is normalized to e ^ 10^°eV. The particles vortons in¬ 
teract with are essentially quarks composing hadrons, so 
the mass m can be taken to be that of the proton (at 
this level of approximation, the mass difference between 
a quark and a proton is negligible). With these num¬ 
bers in mind, Eq. (^ transforms to 7 ^ ^ 10^^ and the 
cosmic string scale d) m ~ 10® GeV which, surprisingly 
enough, turns out to be right at the closure limit given by 
Eq. (Q). It now remains to calculate the effective cross- 
section between a vorton and a quark in order to know 
the expected flux and check that it is compatible with 
the current observational limits. 


II. CROSS-SECTION VORTON-HADRON 

There is no privileged model for cosmic strings, and 
therefore neither is there for describing vortons. Thus, 
we can at best evaluate rough orders of magnitude for 
the cross section we are looking for. There are however 
various levels of approximation at which a vorton can be 


described. As mentioned earlier, the characteristic length 
scale associated with a vorton configuration is expected 
to be a hundred times larger than its thickness, so the 
most obvious description of a vorton is that of a classical 
circle. This is unfortunately of absolutely no use when 
one is looking for the trapped E particles along the loop. 
Hence, as a second level of approximation, we shall con¬ 
sider a vorton to be a torus like configuration in which 
the field S feels a confining potential. This is beyond 
the scope of the present calculation and is left for further 
investigation . So let us describe a vorton as a sphere 
of radius R where in fact the mass of E is a radially 
dependent function. In order to try and recover the cir¬ 
cular geometry, the confined field, expanded as it should 
in the eigenvectors of angular momentum (spherical har¬ 
monics), will be described only by those high values of 
the total angular momentum (to take into account the 
fact that J » I) as well as its projection along a fixed 
axis. 

The basis for the possible bound states of E will there¬ 
fore be given as the set of quantized solutions of the 
Klein-Gordon equation 

(□-HM2)E = 0, (9) 

where the mass M = M{r). The effective radius R is 
adjusted in such a way that the geometrical section of 
the corresponding sphere ttR'^ is that of the equivalent 
ring if it hits a particle face-on, i.e. ttR? = 7 r[(r^, -I- 
— [vy — r^)®], or i?® = 4r^,r^. In these formula, Vy 
stands as before for the vorton radius, while ^ m~^ 
is its thickness, i.e. R ^ Z^/'^mT^. 

In the particular example which was used to calcu¬ 
late numerically the curve on Fig. I, the mass function 
in Eq. © is taken as M{r) = mQ{R — r), and the 
trapped field is assumed the separated form E(a;“) = 
Uneir)YiM{0, , with Yim{S,4>) a spherical har¬ 

monic and the radial function provides the quantized en¬ 
ergy levels for each value of the angular momentum by 
imposing equality of the logarithmic derivative of the nor¬ 
malizable solutions at r = i? (in the simplified model we 
investigate, those are spherical Bessel and Hankel func¬ 
tions of the first kind). As discussed above, we shall 
restrict our attention to AI ~ ~ Z® ~ 10 "* in order to 

take into account, at least partially, the effectively circu¬ 
lar (as opposed to spherical) symmetry. For the actual 
calculations of effective cross-sections, we shall adopt the 
same normalization convention for all the states, i.e. the 
bound and free states (including hadron states) are co- 
variantly normalized at 2E particles in the volume V for 
a state of energy E {uJniM in the case of a bound state). 

We now calculate the cross section ay for an inci¬ 
dent electromagnetically charged particle, a hadron say 
for definiteness, in the frame of the loop, to interact in- 
elastically with a vorton. We consider the situation in 
the rest frame of the loop in which the incident hadron 
with momentum pin = {Ei , pi) hits the E particle bound 
state characterized by its angular quantum number £i 
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and bound energy oji to yield an outgoing hadron with 
momentum pout = iE 2 ,P 2 ) and a final state for the S 
particle which could in principle be either another bound 
state with quantum number if or a free particle state 
with momentum kf = (u! 2 , k); in the notation of the pre¬ 
vious section, one has AE = uj 2 — coi. In both cases, 
we assume the decay rates dTj of the resulting products, 
namely either the excited state or the unstable particle 
itself, to be much smaller than the energies of the parti¬ 
cles so that the cross-section for producing those unstable 
states effectively factorizes as da = dajdTj/Tj. Inte¬ 
grating over the various decay possibilities (f dF^ = F^) 
shows that it suffices to calculate only the intermediate 
cross sections da, and sum them up. Note also that be¬ 
cause our vorton model is very crude, we cannot expect 
more than a rough order of magnitude for the cross sec¬ 
tion, a detailed quantum examination of the same process 
being held for another work . This is the reason why 
we now restrict our attention to the ionisation process, 
including the possibility of final bound state by consid¬ 
ering various energy eigenvalues. 

The interaction is assumed to be only electromagnetic. 
Hence, if ipin and "tpout are the wave functions of the 
hadron under consideration before and after the collision, 
the hadronic current we are interested in reads 

( 10 ) 

a similar definition holding for the S current where 
the “in” state is a bound state, while the “out” state is 
a free S particle. With the electromagnetic field the 
total Hamiltonian for describing the collision is then 

H = + H^+ H^, ( 11 ) 

where the self interaction Hamiltonians and con¬ 
tain in principle all the information about the internal 
structure on the interacting particles. In practice, this 
means that we can work with eigenstates of both these 
self-interacting Hamiltonians to take into account exactly 
the vorton structure (this is supposedly done by calcu¬ 
lating the bound states) as well as the strong interaction 
effects binding the quarks in the hadron together; thus, 
both protons and neutrons can interact with a vorton, a 
point which is understandable by noting that the vorton 
being much smaller than the hadron, it interacts essen¬ 
tially with the quarks. Neglecting these corrections of 
order unity, and working in the Lorentz gauge d,iA>^ = 0, 
the amplitude for this process is easily calculated semi- 
classically as (note we do not consider the fermionic de¬ 
grees of freedom of the hadrons) 

{Pin,ii\Pout,k{if)) = i J d'^xg^‘'Jl^{x)q~'^J^{x) ( 12 ) 


in plane waves and imposing suitable commutation rela¬ 
tions (anticommutation to describe properly the hadron 
states) iQ. We recall that our purpose here is simply 
to derive a rough order of magnitude estimate for the 
interaction cross section, which explains why we perform 
only a semiclassical analysis and neglect the fermionic be¬ 
haviour of hadrons. Use of Eq. (12) shall be made later 
for deriving the details of our oversimplified spherical 
model, but for the time being, let us derive the expected 
order of magnitude of the interaction cross-section. 

The total cross-section, being a sum over all the con¬ 
fined particles of terms like Eq. & squared is seen to 
have a factor Z^e'^, with e the electromagnetic coupling 
constant. Then, on dimensional grounds, it can be noth¬ 
ing but proportional to {AE)~‘^. Thus, one ends up with 


2„4 


Z^e 


= 


{AE} 


~ 10 


-32 


10^ GeV 


m 


7 ^ 

^100 


cm 


(13) 


where ^loo = 2^/100 and the function is a dimension¬ 
less quantity thats need a specific modef to be evaluated. 
In the spherical symmetry approximation and for a few 
energy levels, it is exemplified on Fig. ^ The line and 
continuum spectrum characteristic of bound state inter¬ 
actions is clearly shown on the figure, as well as the exis¬ 
tence of an energy threshold below which no interaction 
takes place. 



FIG. 1. Expected qualitative form of the cosmic ray spec¬ 
trum for vortons: it consists of a line spectrum followed by a 
continuum. This function J^(Ai?/m) [see Eq. (@)] modulates 
the cross-section and was calculated using 6 energy levels with 
I = 10“^, taking everything numerically into account. 


where q = Pout —Pin is the exchanged momentum. At this 
point, a completely quantum field analysis could be made 
by expanding the E operator in bound and free particle 
creation and annihilation operators and the hadron field 


To conclude this section, we consider the probability 
that a vorton interacts in the way discussed above in the 
atmosphere. Using the values derived in the previous sec¬ 
tion for the mass and charge of the vorton shows that the 
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characteristic ionisation cross section is typical of neu¬ 
trino interaction at these energies ||l^ . The probability 
we are looking for is therefore, using pAtm ~ 10 ^ g-cm”^ 
for the mean atmospheric depth 

ay = UyPAtTnlmp ~ 10“^ - 10“^Zfoo. (14) 

This quantity we shall use later to compute the expected 
flux of this type of events on earth. 

III. ACCELERATION MECHANISM FOR 
VORTONS 

In this section and the following, we turn back to or¬ 
dinary units in which c and h have their usual values. 

Accelerating a particle of charge Ze like a vorton to 
energies larger than 10^° eV by means of an electric 
mechanism requires a potential difference larger than 
10^°/Z Volts. There are basically two classes of as- 
trophysical objects in which such potential differences 
may be found, namely pulsars and accreting black holes 
(BH). In fact, potential differences as high as 10^® Volts 
are known to exist in the magnetosphere of young pul¬ 
sars and in that of the accretion disks around massive 
(10^ - 10 ®Mq) Kerr BH. 

Pulsars as sources of high energy vortons can be im¬ 
mediately ruled out since vortons are not expected to be 
present at the surface of neutron stars. Their mass is so 
high, as we have just seen, that they would sink toward 
the center of the neutron star. We shall therefore for now 
on consider the case of accreting BH, supposing those to 
be the power engine for active galactic nuclei (AGN) ra¬ 
dio galaxies and quasars (QSO). The underlying idea is 
the following: radio jets, X and 7 rays emissions are due 
to the acceleration of electrons, positrons and protons by 
electrostatic fields generated by a Blandford type mecha¬ 
nism near the horizon of the BH. The model is based 
on the assumption that a fraction a of the total mat¬ 
ter is made of vortons jl^, which behaves essentially as 
cold dark matter. Therefore, one can expect their spa¬ 
tial distribution to look like that of ordinary matter. The 
main idea is then that in an accelerating object, what¬ 
ever it is, the same fraction a of vortons is accelerated by 
the same mechanism that works for the protons. How¬ 
ever, in reasonable models, the total mass of the vorton 
is My Zm where Z is the charge carried by the vor¬ 
ton and m the scale of symmetry breaking at which the 
strings were first formed. As shown in Sec. a Lorentz 
acceleration factor of 10 ® — 10 ® is actually sufficient to 
reproduce the data. This means that a very basic accel¬ 
eration by means of an electric field is enough since at 
these energies, losses in radiation can be considered neg¬ 
ligible. In what follows, we shall justify the above model 
to show that it encompasses no major difficulties. 

It should first be clear, and we want to emphasize that 
point again, that no exotic mechanism is required for 
the acceleration (even though the accelerated particles 


themselves may be considered exotic). Indeed, the accel¬ 
eration mechanism that produces ultra-relativistic jets in 
radio galaxies and X and 7 rays in AGNs and QSOs is 
enough to accelerate vortons to energies up to 10^® eV 
and higher. We want to point out that electrostatic helds 
are very likely to be responsible for the acceleration of 
particles in jets of radio galaxies and for the generation 
of high energy photons (with energies exceeding I TeV). 
Although the origin of such electrostatic helds is not com¬ 
pletely understood yet, they may well originate through 
the quite appealing Blandford-Znajecs mechanism [ p4| . 
Let us hrst describe shortly such a mechanism. 



FIG. 2. The electromagnetic field configuration surround¬ 
ing a Kerr BH in the poloidal Blandford Znajecs mechanism. 
The circle represents the event horizon of the BH (r = r+) 
and the shaded area stands for the accretion disk. 

Consider a magnetized rotating neutron star. Its sur¬ 
face is supposed to be a perfect conductor, so that in the 
rest frame of the surface of the neutron star, the tangen¬ 
tial component Eg of the electrostatic held must vanish, a 
condition which, when written in the inertial frame takes 
the form 

= (15) 

where n, B and i?* are respectively the angular velocity, 
surface magnetic held and radius of the neutron star. 
This gives very energetic electric held lines along which 
charged particles can be easily accelerated. 

An analogous mechanism for generating electric held 
lines works for Kerr BHs, with the unfortunate differ¬ 
ence that a bare BH cannot have a magnetic held, the 
latter being supplied by the accretion disk via a dynamo 
mechanism. The surface of the BH behaves like a ro¬ 
tating conductor with angular velocity H = dc/M, with 
a and M the Kerr solution parameters somehow identih- 
able with the total mass (M) and the angular momentum 
(for 0 < d < I) of the source. Now because of the bound¬ 
ary conditions on the surface of the BH, an electrostatic 
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field is created which yields a potential difference given 

by [H 


Ay = 


aB 

H 




(16) 


where B is the poloidal component of the magnetic field 
generated by the accretion disk and D the typical length 
scale of the electrostatic field, D ~ M, i.e., Ay ~ aBM 
(see Fig. |). 

The magnetic field is estimated to be 10^ gauss for the 
most active galaxies. In fact, under the hypothesis that 
the magnetic energy density is in equipartition equilib¬ 
rium with the thermal energy of the disk, one finds 

/ r \ 1/2 2 

B = j _ ~ 4 X gauss, (17) 


main mechanism allowing acceleration of protons at en¬ 
ergies ~ 10^° eV, while protons in the near zone of the BH 
cannot be accelerated at so high energies because of the 
losses due to the unfavorable environment (synchrotron 
and Compton losses) . Vortons on the other hand do 
not suffer from this drawback since because of their high 
mass, they have a Lorentz factor which is at most of the 
order of 10 ® as we have seen above. 

The synchrotron energy loss per unit length is 

^ = (22) 

ax o p 

with p the curvature radius of the vorton trajectory 
along the magnetic field lines. If we take p equal to the 
Schwarzchild radius = 2GM/c^ of the BH, we obtain 
the total energy loss 5E as 


where is the Eddington luminosity 


^ AnGcMnir, -i 

^ ~ 10 ^®Mg erg • s \ 


(18) 


TOp being the mass of the proton, CTj, the Thompson cross- 
section 1^ and Mg = (M/10®M©) being the BH mass in 
units of 10® solar masses. Taking into account Eq. (16), 
we then end up with a potential difference given by 


Ay ~ 1.7 X 10^°aMi/^ Volts 


(19) 


or, equivalently 


Ay = \j —a ~ 5 X 10 ^®aL 45 Volts, 
c 


( 20 ) 


L 45 being the luminosity in units of 10^® erg/sec [ p5[ . 
Estimates of B [Eq. (^] and Ay [Eq. (^)] are very 
rough and depend crucially on the conversion efficiency 
of gravitational energy of the disk into electromagnetic 
energy. As already stated, the parameter a is less than 
unity and is related with the angular momentum J of the 
BH through 


J=-GM^, (21) 

c 

where M is now expressed in ordinary units. Eor old 
AGNs and QSOs, a can be larger than 0.9. 

So high a potential difference is very efficient to accel¬ 
erate particles and can thus be the source for the elec¬ 
tromagnetic energy of AGNs, QSOs and radio galaxies. 
Similar powerful sources of electromagnetic energy can 
be found in the magnetosphere of the accretion disks and 
by radiation accelerated winds from the disk. Interested 
readers will find details and useful references in the re¬ 
view article 

Hot spots in radio galaxy jets are also good candi¬ 
dates as acceleration regions of very high energy parti¬ 
cles. Shock waves in the jets of radio galaxies are the 


“ = 3Glf‘ ■*'' 

= 5.2 X 10"®®ZMM8"S'^ erg 


(23) 


where A is the acceleration path expressed in 
Schwarzchild units. For Z = 100 and 7 = 10®, the syn¬ 
chrotron losses are ~ 10 “"^ erg (with A= 1 ), a very small 
quantity indeed with respect to the required final energy 
of 10 ® erg. 

A similar calculation can be performed for the Comp¬ 
ton losses; let us first compute the Thompson cross sec¬ 
tion cr^'* of a vorton: 


a 


Th 

V 


Stt 


{Zef 

ZMIA 


= 1.6 X 



2 

cm^ 


(24) 


(recall My = Zm). It is now easy to compute the mean 
free path of a vorton in the photon bath around a BH. 
For a given luminosity L, the photon density nph reads 

where R is the radius of the source and is a typical 
frequency of the radiated spectrum. Again fixing R to 
be the Schwarzchild radius of the BH, we have 

Uph = 2 X 10 ^®L 45 Mg~^ ( ) pbotons/cm^, (26) 

V ^ph ) 

so that the vorton mean free path = {nph<Ty'')~^ is 

(27) 


which is much greater than the typical acceleration length 
4 ~ 3 X 101®M8 cm. 
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In brief, the acceleration of vortons in the hot spots 
of radio galaxies is more favorable than the same mecha¬ 
nism applied to protons because synchrotron and Comp¬ 
ton losses scale respectively like 

(^) “ (synchrotron) 


and 


(Compton). 

In what follows, we postulate that vortons that are uni¬ 
formly distributed in the accretion disk around a BH can 
reach the region of high electrostatic field. Taking into 
account Eq. ( 0 ), we see that they can gain energies of 

= m‘^°ZaL 45 eV, (28) 

so that even active galaxies with luminosities as low as 
1043 ej-g/g pan accelerate vorton to 10^° eV if the angular 
momentum of the BH is large enough that a ^ 1. 

As already stated, acceleration at these energies of 
vortons can happen also in the radio galaxy hot spots. 
The possibility of accelerating them near the central BH 
means that under the hypothesis that monopolar electro¬ 
static and rotating BH is the power engine of all active 
galaxies [^, Seyfert I, H, QSO and AGN are potential 
vorton acceleration sites, thereby tremendously increas¬ 
ing the number of sites as compared to the proton case 
since there is no GZK cutoff in their case, as we shall now 
see. 


IV. ENERGY CONSIDERATIONS AND 
PROPAGATION. 

The flux of cosmic rays with energy higher than 10^° 
eV is estimated of the order of 10“^° cm“^-s“^-ster“^, 
i.e. a flux 10 “^® cm“®-s“^, so the energy flux of particles 
having energies E = E 20 x 10®®eV is 

= 10 ~^^E 2 o erg • cm“^ • s“^. (29) 

If Vy is the probability that a vorton interacts with the 
atmosphere and if all the particles with energies greater 
than 10®° eV are vortons, then the actual flux (as opposed 
to the observed one) is 

$ = 10 -^E 2 oV 5 erg • cm”® • s-\ (30) 

with V 5 = Vy/10~^ (see Sec. i). In order to achieve 
this flux, extragalactic sources up to a distance D must 
supply a power W 


W = 10‘^^E2oV5 


/lOO Mpc 

1 ^ 


2 

erg • s 


(31) 


Assuming 10^ potential sources (Sy I, I, QSO) means 
that each source must supply a power 

W = 10^^E2or5 erg • s'^, (32) 


power too uncomfortably close to the averaged power ra¬ 
diated by active galaxies. Note the neutrino hypothesis 
suffers even more of the same drawback. 

Until now, we have assumed the typical mean free path 
of a 10®° eV vorton to be much larger than the Hubble 
radius, with the meaning that a vorton can reach us, 
even if emitted at large redshifts. We shall now prove 
that point by estimating the energy losses the vorton has 
to experience on his way. 

We start with synchrotron radiation due to intergalac- 
tic magnetic field B since this is supposed to be the most 
efficient mechanism. The energy loss per unit time is 


dt 3 A2 ^ ’ 


(33) 


in a mean magnetic field B = if the Larmor 

radius of the vorton is less than the correlation length of 
the magnetic field B. The typical slow down timescale 
for the vorton is thus given by 


, 3 My 
^ 2E{B^){Ze)^ 


= 4.6 X 10 


59 


B 


E. 


-1 (. 


m 


10 “®gauss J 20 lO^TO 


(34) 


which is clearly sufficiently large that there is not even a 
cosmological cutoff for vortons. 

A similar conclusion holds for the energy losses due to 
the microwave background. Note however the essential 
difference in this case between protons and vortons: if 
protons were pointlike particles with no internal quark 
structure, they would be able to propagate almost freely 
in the microwave background, even at these energies, and 
they decay dominantly because the binding energy be¬ 
tween quarks and gluons is much less than the proton 
mass. Vortons, on the other hand, not only have en¬ 
ergy levels which are much higher than protons, but also 
they propagate with a relatively low velocity (7 ~ 10 ®). 
Therefore, vortons, once accelerated can arrive on earth 
with undegraded energy: supposing they were acceler¬ 
ated at the time of formation of the galaxies {z ^ 1 ) [^ , 
they would by now have lost a mere factor of two because 
of cosmological redshift. The interesting point in that 
observation is that for z ^ 1 , the ratio between active 
and normal galaxies is about 0.1 |^ , ^ . Keeping that in 
mind, we can now calculate the high energy vorton den¬ 
sity under the hypothesis that vortons were accelerated 
at the birth of galaxies. 

If one assumes that all the unexplained cosmic rays 
with energy larger that 10®®eV are made of vortons, their 
density will be 
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= 4.2 X V • cm-\ (35) 

C I y 

where (and $20 = cm“^.s“^.ster“^) is the 

high energy cosmic ray flux. 

The total number of vortons Ny in a sphere of radius 
corresponding to 2 = 1 is 

Ny=Sx (36) 


and their energy is 10^°(1 + z)Ny eV = 6 x 

10 ®®/i“^$ 2 o/ 7^5 ergs, where h is the Hubble constant Hq 
in units of 75 km-s^^-Mpc”^. Such an amount of energy 
can be released by a number of active galaxies 


Na = 


2 X 10®. 


_3 4'20„_i 


V 5 


Tg T43, 


(37) 


with Tg the duty time of the galaxies in millions of years 
and £ the ratio between the power used to accelerate vor¬ 
tons and the electromagnetic luminosity; note this latter 
parameter is not constrained and can in fact exceed unity. 

The baryonic mass in the sphere with z = 1 is = 
5.8 X g, so that the total number of galaxies is 

fVc - 3 X (38) 

where Mm is the galaxy mass in 10^° solar masses. From 
Refs. we know the ratio between active and nor¬ 

mal galaxies to be of order 10%, so the required efficiency 
is, expressed as a function of solar masses accreted by the 
BH per year M, 


£ 


Fs 



(39) 


Finally, let us compute the required number of vortons 
present in the universe in the form of dark matter. Let a 
be the ratio between the total vorton and baryon masses 
in the universe. The ratio r] between vorton and baryon 
number densities is 

and assuming the baryon density to be 10% of the critical 
density {rn, = 6.6 x“^ h? baryons/cm®) we finally obtain 


a = 6.4 X lO"'^ 


‘h20 


^100 



(41) 


This result is understandable in two ways. Either all 
vortons present in AGNs are indeed accelerated, in which 
case vortons are, according to this calculation, expected 
to represent roughly less than a thousandth of the matter 
in the universe, or, conversely, vortons do fill the universe 
so that the actual value of a should exceed unity. In 
the latter case, our calculation reveals that either the 
interaction probability is lower than what we evaluated 
above, or only a small fraction of vortons get accelerated. 
In any case, it should be clear that our model provides a 
very high energy cosmic ray flux which can be made to 
agree with observations. 


CONCLUSIONS 

We have exhibited a model for explaining extremely 
high energy cosmic rays that have recently been observed. 
More of these events are expected to be observed in the 
near future by the Auger Observatory |l^]. We propose 
that they are bound states of very massive particles in 
vortons, i.e. loops of superconducting cosmic strings sta¬ 
bilized by a current, that are are freed by inelastic col¬ 
lision with atmospheric hadrons. The model, in its sim¬ 
plest form, has only one free parameter, namely the en¬ 
ergy scale at which the current sets up in cosmic strings. 
This could easily be extended to take into account a pos¬ 
sible difference between this mass scale and the string 
energy scale itself . In this particular framework how¬ 
ever, the mass scale is constrained by requiring that vor¬ 
tons do not overfill the universe, a bound that is almost 
saturated by the demand that high energy cosmic rays 
are made of vortons. This interesting coincidence could 
imply, if verified, that a non negligible fraction of the 
dark matter in the universe consists of vortons. 

Let us here summarize our findings concerning the 
model. First, as mentioned earlier, it essentially has only 
one free parameter, namely the mass scale at which the 
superconducting current sets up in the core of the cosmic 
strings. As it turns out, demanding these objects to be 
candidates for the few 10®® eV events through interaction 
with atmospheric protons completely fixes this parame¬ 
ter to 10® GeV. There doesn’t seem to be any way out of 
this prediction, which renders the model very falsifiable. 

Another point worth mentioning is that the vorton- 
proton (or neutron) cross section is quite weak, giving, 
at these energies, a total interaction probability between 
10“® and 10“^ with the earth atmosphere. Note again 
that this probability depends on nothing but the fixed 
energy scale, and that, in order to fit the observed data, 
it implies a high energy vorton flux whose numerical value 
is also therefore determined. 

No high energy cosmic ray model is complete without 
specifying acceleration and propagation processes. In our 
case, acceleration is performed very simply by kicking the 
vortons (at least those which have been ionized somehow 
beforehand) with the high electrostatic fields that are 
expected to be present in AGNs. The required quantity 
of energy turns out to exist in these objects and because 
of the relatively low amount of acceleration required (7 < 
10®), losses are negligible and the vortons can escape the 
acceleration zone. This is, to the best of our knowledge, 
the most efficient mean of extracting 10®® eV in a single 
particle out of any astrophysical object. 

Similarly to the fact that energy losses are negligible 
in the region of acceleration, and for the same reason, 
the propagation in intergalactic medium is done almost 
without any collision, thus making this medium effec¬ 
tively transparent to very high energy vortons. As a re¬ 
sult, there is no reason for a cosmological (GZK) cutoff 
and vortons can come from as far as a redshift of a few. 




This is satisfactory because the ratio of active to normal 
galaxies increases as one goes farther away, and the final 
vorton flux we end up with thanks to this fact is actually 
quite large and indeed requires only a very small fraction 
of the total mass density in the universe to consist in vor- 
tons (of the order of 10“'^) if all vortons are ionized (and 
therefore accelerated), or, if one believes the standard 
vorton predictions [Q that for a mass scale of 10® GeV 
they should be the dominant part of the dark matter, 
then one just needs the same tiny fraction of vorton to 
be accelerated. The absence of any cutoff also means we 
predict a spatial isotropy together with a correlation with 
active galaxies. 

Interacting with the atmosphere with a very low prob¬ 
ability (and then presenting a threshold followed by a 
line spectrum, presumably unobservable unfortunately), 
vortons are here predicted to give rise to mostly hori¬ 
zontal air showers, just like the neutrinos. Let us stress 
however that until now, it has not been possible to find 
any astrophysical mechanism that would give neutrinos 
such high energies. Hence, although our vortons are in¬ 
deed hypothetical in the sense that we don’t know what 
is the actual theory that describes physics beyond the 
electroweak scale at which we expect to find them, it 
should however be clear that they do not, apart in their 
existence itself, imply any new mechanism of any kind. 

Finally it is worth pointing out that the acceleration 
mechanism we have proposed is just one amongst a few 
other acceleration mechanisms already proposed in the 
same context such as the Fermi mechanism. Vortons can 
be accelerated as well across the shocks of the jets of the 
radio galaxies. We prefer however the direct acceleration 
mechanism because of its high efficiency and therefore 
the energy budget constraints can be more easily fulfilled. 
Note that this is not the case for the neutrinos: in fact, as 
already stated, 10^®eV neutrinos have interaction cross 
sections (with hadrons) that are roughly equivalent to 
vorton’s. Consequently, it is highly questionable whether 
the energetic constraints required to accelerate the parent 
charged particle of the neutrinos in the lobes of radio 
galaxies can be satisfied. 

To end up this conclusion, let us compare our expecta¬ 
tions with those of the Auger Observatory ||^, a project 
specifically designed to accumulate more statistics on cos¬ 
mic rays with energies in excess of 10^® eV. The Auger 
Observatory will have an energy resolution less than 20 % 
which is presumably going to be insufficient to actually 
allow definite conclusion about the line features. In spite 
of this difficulty, it is however not impossible that part of 
the line signal might be found using correlation function 
techniques after a power-law substraction will have been 
made in the data: no significant feature should be left by 
almost all the competing candidates. Another character¬ 
istic of the Auger Observatory concerns its ability to iden¬ 
tify the primary of an air shower through accurate mea¬ 
surement of the shower maximum (AX^ax ^ 20g-cm“®). 
In a bound state model such as the one we propose here, 
the line part of the spectrum should be initiated by ra¬ 


diative decay of the excited vortons, hence the primary 
should be a photon. Then once the continuum is reached, 
some other particle (the E in our model) plays the part 
of the primary, initiating a shower whose maximum is 
expected somewhere else. It is not clear yet, because 
very model dependent, whether the Auger Observatory 
accuracy in this measurement will be sufficient, but it 
is also a firm prediction of any bound state model. Fi¬ 
nally, according to our calculation, the interaction prob¬ 
ability with the atmosphere is very low, with the result 
that many horizontal showers are expected. Contrary to 
most other cosmic ray detectors, the Auger Observatory 
will be very efficient to detect those, being also under¬ 
standable as a neutrino detector with an acceptance of 
10 km^-sr water equivalent for horizontal air showers at 
10^®eV 1^^. Finally, although the angular resolution of 
the two detectors (giant array and air fluorescence) in the 
Auger Observatory should be of the order of less than 
2 °, some events will be observed in hybrid mode by both 
with a resolution of 0.3°, a precision that is expected to 
be enough to conclude on the isotropy of the cosmic ray 
sources. 

All these facts lead to the conclusion that our model 
has a very high potential for being either confirmed or 
ruled out shortly after the Auger Observatory will be 
started. 
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